Phase transition in the collective migration of tissue cells: experiment and model 
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We have recorded the swarming-like collective migration of a large number of keratocytes (tissue cells ob- 
tained from the scales of goldfish) using long-term videomicroscopy. By increasing the overall density of the 
migrating cells, we have been able to demonstrate experimentally a kinetic phase transition from a disordered 
into an ordered state. Near the critical density a complex picture emerges with interacting clusters of cells 
moving in groups. Motivated by these experiments we have constructed a flocking model that exhibits a contin- 
uous transition to the ordered phase, while assuming only short-range interactions and no explicit information 
about the knowledge of the directions of motion of neighbors. Placing cells in microfabricated arenas we found 
spectacular whirling behavior which we could also reproduce in simulations. 

PACS numbers: 87.17.Aa, 87.17.Jj, 87.64.Rr 



I. INTRODUCTION 

The collective motion of organisms is a spectacular phe- 
nomenon sometimes involving huge schools of fish, thou- 
sands of birds exhibiting complex aerial displays 111], herds of 
quadrupeds and even bacteria producing fractal colonies Q2Q 
or amoeba assembling into rotating aggregates {H. Recently, 
sperm cells were also demonstrated to form self-organized 
vortices [4J]. In addition to being a common mechanism by 
which organisms self-organize, a deeper understanding of the 
simultaneous adjustment of the velocities of many moving ob- 
jects has important potential applications ranging from the 
swarming of distributed robots exploring new territories [5] 
to the healing of wounds related to the coherent migration of 
epithelial cells jg]. 

Although widely observed in nature, collective motion is 
less accessible for experimental investigations under labora- 
tory conditions. Becco et al. |7] have presented interesting 
but only qualitative results. A well controlled series of experi- 
ments aimed at characterizing the nature of the transition from 
a disordered to an ordered phase in velocity space is so far 
lacking, even though such experiments would be useful both 
in providing a quantitative reference for further experimen- 
tal studies and in prompting more realistic models for group 
behavior. While laboratory observations have been scarce, a 
number of models based on self-propelled particles have been 
developed recently (see, e.g., Refs. Jill [Sill [II HQ) 
to describe collective motion. In broad terms these mod- 
els fall into two categories, those which describe the onset 
of collective motion as a transition to an ordered state in 
a large noisy system with simple interactions between the 
particles [8l |s£ llfJl . and those which employ more compli- 
cated - consequently more realistic - interactions and focus 
on the emer gent collective dynamics of groups of finite size 
ITll [l2l [l3l Il4ll . A model, however, that combines simple 
short-ranged interactions, detailed realistic dynamics and a 
well defined kinetic transition in a large noisy system has not 
been described so far. 

In this paper we present experimental results concerning the 
collective migration of a large number of tissue cells (kerato- 



cytes) using long-term videomicroscopy. As the overall den- 
sity of the migrating cells is increased, we observe a kinetic 
phase transition from a disordered (low density) state into an 
ordered (high density) state, in which most of the cells move 
in a direction approximately agreeing with their average di- 
rection of motion. Just below the transition a complex picture 
emerges with interacting groups of cells moving in random di- 
rections. Motivated by these experimental results we develop 
a flocking model, which, in contrast to previous models, con- 
siders a minimal realistic interaction that assumes no explicit 
averaging of the directions of motion, while also exhibiting 
a transition to the ordered phase. Numerical studies indicate 
that this transition is continuous and belongs to the same uni- 
versality class as the model of Vicsek et al. HI]. These results 
compel us to imply that the experimental transition described 
is continuous as well. 



n. EXPERIMENTAL SETUP AND RESULTS 

Our experimental setup, consisting of a home developed 
fully computer-controlled time-lapse microscope [15] and a 
custom made room temperature CO2 mini-incubator allowed 
us to carry out long-term videomicroscopy of keratocytes to- 
gether with a quantitative analysis of their motion. We col- 
lected 2-4 fish scales from living goldfish {Carassius auratus) 
with tweezers, and placed them external side up in a 35-mm 
circular Petri dish similarly to Ref. 11611 . Scales were kept 
in the incubator overnight to allow epidermal keratocyte cells 
to migrate out from the scales. Before time-lapse microscopy 
scales were removed, and the cells remaining in the Petri dish 
were treated with phosphate saline buffer and/or trypsin to ob- 
tain cultures of varying density (both compounds reversibly 
weaken cell-cell and cell-substrate connections, allowing the 
removal of a controllable percentage of the cells). During the 
subsequent typically 24-h long time-lapse microscopy exper- 
iments we monitored the motion of live keratocyte cells (tak- 
ing pictures in several fields of view at frequency of 1 shot 
per minute) in cultures of varying density ranging from sparse 
ones with very low cell densities to confluent ones with nearly 
complete coverage. We observed a relatively sharp transition 
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FIG. 1: Phase contrast images showing the typical behavior of cells for three different densities, (a): 1.8, (b): 5.3, (c): 14.7 cells/lOOxlOO 
frni 2 . We observed that as cell density increases cell motility undergoes collective ordering. The speed of single cells is higher than that of 
cells moving in coherent groups. Scale bar: 200 /im. (d)-(f): Velocity of cells. Scale bar: 50 ^im/min. Online supplemental material fnll 
contains corresponding videos. 
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FIG. 2: Order parameter V is shown as a function of normalized 
cell density. Cell density was normalized with the maximal observed 
density of 2.5 x 10~ 3 cells//im 2 and error bars indicate the standard 
error of the density and order parameter. 



from random motility to an ordered collective migration of 
dense islands of cells as the density was increased. Fig. Q] 
shows the typical behavior of cells for three different densi- 
ties. The technical details of the experiment are described in 
the online supplemental material fllTfl . 

In order to quantify the level of coherence of the migration 
pattern of cells -the individual motion of which are usually 



described as a persistent random walk HHl - we calculated the 
velocity of 20-30 cells in each experiment from their displace- 
ment between frames. This process yielded the Vi(tfc) veloc- 
ity vectors of cells, where i is the index of the cell, tk — fcAi 
is the time elapsed from the start of the cell's trajectory and 
At denotes the time difference between frames. We define as 
the measure of coherence of motion an order parameter equal 
to the time average of the sum of the normalized velocities 
divided by the number of cells measured. Thus, the order pa- 
rameter is 



V= ^ 



N 

£■ 

i=l 
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where N is the number of evaluated cells. Fig. [2] displays the 
order parameter as a function of cell density. Our measure- 
ments were carried out after the cells had had sufficient time 
to migrate out of the scales and achieve a quasi-stationary mi- 
gration. V was calculated by averaging over the observation 
time and its standard error was calculated by dividing the stan- 
dard deviation by y/Al — 1, where M is the number of ana- 
lyzed snapshots. We attempted to survey the extent of finite 
size effects, by calculating the order parameter and its error, 
while considering only half the number of available cells. We 
found no significant deviation. Cell density was measured lo- 
cally, i.e. in the field of view, every 30 minutes. Both its mean 
value and the standard error of the mean is presented in Fig. 2. 
A sharp increase can be observed around 5 x 10~ 4 cells//im 2 . 

We conclude that a kinetic phase transition takes place from 
a disordered into an ordered state as cell density exceeds a 
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FIG. 3: Computer simulations were performed of the system described in the text for different densities. These simulations showed a transition 
to the ordered phase similar to that seen in experiments, (a) Typical behavior of cells is shown for three different values of the normalized 
number density p — p/p ma x, with p max ~ 2, which is approximately the density where gaps disappear, and the cells reach tight packing in 
simulations, (b) As cells moving in different directions come into contact adhesive inter-cellular forces act to align ru and rij, resulting in an 
effective averaging of self-driving directions. For videos of simulation runs see the online supplemental material H 1 711 - 



relatively well-defined critical value. Our experiments suggest 
that short-range attractive-repulsive intercellular forces alone 
are sufficient to organize motile keratocyte cells into coherent 
groups. 



III. MODEL DESCRIPTION AND RESULTS 

To interpret the above phase transition-like ordering phe- 
nomenon we constructed a model that takes into account the 
specific features of the experimental system and is able to re- 
produce experiential behavior. In our model individual model 
cells (self-propelled particles) move forward in a well defined 
direction with constant speeds. The noisy nature of the pro- 
cesses which generate cell locomotion is taken into account 
by considering the direction of self-propulsion of model cells 
to be noisy. Intercellular forces through which model cells 
interact are considered to be short-ranged, as between kera- 
tocytes they are the result of direct physical contact. Further, 
regarding interactions between keratocyte cells, it is obvious 
that explicit averaging of the directions of motion employed 
in previous models is not realistic. Tissue cells forming co- 
herently migrating groups are unable to explicitly adjust their 
direction of motion to the average velocity of their neighbors, 
collective motion must emerge solely as a result of direct cell 
to cell interactions (forces). To model the emergence of col- 
lective motion without such explicit averaging, we consider 
self-propelled particles (model cells) that attempt to adjust 
their direction of motion toward the direction of the net-force 
acting on them. 

The 2-dimensional motion of model cell i £ {1, N} with 
position Yi{t) is described by the overdamped dynamics: 



dr<(t) 
dt 



N 



uoni(t) +M^ F ( r ii r j)- 



(2) 



propelling velocity of magnitude vq in the direction of the unit 
vector rij(t) and experiences intercellular forces F(rj,rj). 
The direction of the self-propelling velocity rij(<), described 
by the angle 9f(t), attempts to relax to Vj(t) = dr.;(t)/di 
with a relaxation time t, while also experiencing angular 
noise £: 



d#f(*) 1 . ,( ,s 



(3) 



where e z is a unit vector orthogonal to the plane of mo- 
tion (see Fig. 03) and £ is a delta correlated Gaussian white 
noise term with zero mean, ie. (£(t)) = and = 
ry 2 /12 6(t,t'). We consider pairwise inter-cellular forces 
whose magnitude is a function of the distance between the 
two centers of mass only. Aiming for simplicity we present 
results obtained by using a piecewise linear force function, 
the existence of the transition we describe, however, does 
not depend on the specific form of the function employed. 
The piecewise linear force function we considered was repul- 
sive for distances smaller then R cq ., attractive for distances 
Rcq. < dij < Rq and zero for cells farther apart, i.e. 
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is the value of the 
is the maxi- 
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Thus, each cell with mobility /i attempts to maintain a self- 



-y ~ \n-Tj |' — r» 
maximum repulsive force at dij 
mum attractive force (resulting from adhesive interactions be- 
tween cells) at dij = Rq (see Fig. 03). The values vq = 1, 
fi = 1, t = 1, Rq = 1, R ccl . = 5/6, F adh . = 0.75 and 
■Prop. = 30 were used, the results obtained, however, were 
not sensitive to the particular choice for the parameter values. 
The parameters of the function (@}, the slope of the two linear 
segments (F rcp and i^adh.) an d the equilibrium distance R cq ., 
were adjusted while observing the simulations and comparing 
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FIG. 4: The average value and standard error of V is shown as a 
function of the normalized number density (a) p = p/p max and (b) 
noise rj. Each data point was obtained from at least 10 independent 
simulation runs with N = 1000, with (a) rj — 0.6, p ma x ~ 2 and 
(b) p = 0.6. The insets shows the dependence of InV on (a) (p — 
Pc(t]))/Pc and (b) (r) c (p) ~ ? ))/ J 7c, the slope of the fitted lines can 
be associated with the critical exponents 5 and j3 . The large scaling 
regimes and the similarity of the numerical values obtained for the 
exponents with those found for other models indicate the existence 
of a continuous phase transition in both p and r\. 



them with experimental videos (cf. online supplemental ma- 
terial), in order to achevie best possible agreement with the 
experiments. The two main criteria were: i.) the reproduc- 
tion of the 2-dimensional sheet-like motion, wherein individ- 
ual cells only interact with their nearest neighbors - i.e. cells 
have a well defined volume and ii) the ability of cells to "break 
free" from each other. The first was readily achievable by set- 
ting the slope of the repulsive force segment to be sufficiently 
larger than the that of the attractive segment and the value of 
i? cq . sufficiently larger than Rq, while the second necessitated 
that the maximal value of the attractive segment i^dh. be set 
smaller then vq / p. 

We carried out simulations (Fig. [3^) by solving the system 
of 3iV differential equations (O and ® with periodic bound- 
ary conditions in systems of size Lx L using a fixed time step 
of At — 0.05 Rq/vq, The angular noise term £ was mod- 
eled by choosing its magnitude uniformly from the interval 
[—T]/(2y/At),r]/(2y/At)]. In good agreement with experi- 



ments, running our model resulted in a continuous transition 
to the ordered phase (where all cells move in a common di- 
rection) as the number density of cells p — N/(L/Rq) 2 was 
increased for a fixed value of the noise r\ = 0.6 (Fig. [4J. 
While varying the value of r\ was not experimentally feasi- 
ble, simulations show that a continuous transition to the or- 
dered phase also takes place as r\ is decreased while p is held 
fixed (data not shown, see online supplemental material Il7ll ). 
In both cases we calculated the order parameter ([TJ using the 
velocities Vj(tfc) obtained from the numerical solution of (O 
and ©. V and its standard error were calculated by aver- 
aging over at least 10 independent simulation runs each at 
least an order of magnitude longer than the typical velocity 
autocorrelation time (which was found to be in the range of 
10 - 1000 Rq/vo). The values p c = 0.18 (p c = 0.09) and 
r\ c = 1.18 were found for the critical density and noise from 
simulations of systems with TV = 250, 500, 1000 and 2000 at 
fixed values of rj = 0.6 and p = 0.3 respectively. Based on the 
assumption that our model exhibits a kinetic phase transition 
in the thermodynamic limit analogous to the continuous phase 
transition in equilibrium systems, we proceeded to study it's 
critical behavior, that is, 



V oc [r] c (p) - T]\ 



and 



^cx [p~pM] 6 



(5) 



where (3 and S are the critical exponents and r/ c (p) and p c {rf) 
are the critical noise and density. Analysis of the data yielded 
the values (3 = 0.44 ± 0.08 and 8 = 0.38 ± 0.07, strongly 
suggesting that our model belongs to the same universality 
class as the angular noise model of Vicsek et al. 081], for which 
j3 = 0.45 ± 0.07 and S = 0.35 ± 0.06. 

When two cells approach each other close enough, equa- 
tion (O leads to a gradual alignment of their direction of mo- 
tion. Our model is in this sense similar to other models of 
systems of self-driven particles exhibiting emergent collective 
motion - i.e. flocking -, with the very important difference 
that the particles - the cells - do not directly use informa- 
tion on the movement of others around them to determine 
their own movement. Also, while several other models in- 
clude self-propelling particles, which interact through various 
forces (and the seminal work of Shimoyama et al. ifTTll has 
alignment dynamics similar to ours), they all rely on either 
long-range forces lfTTl[T4Tl . explicit averaging [8] or both lfl3ll . 
The model presented above, on the other hand, combines an 
experimentally motivated minimal dynamics with short-range 
adhesive and repulsive forces, while also displaying a contin- 
uous kinetic transition to the ordered state. 



IV. EFFECTS OF BOUNDARY CONDITIONS 

To investigate the effects of boundary conditions on col- 
lective cell motility, we used square as well as more com- 
plex shaped microfabricated arenas, which kept cells in a well 
defined area. Microstructures were fabricated by UV-curing 
Norland optical adhesive (NOA63) on the surface of glass 
cover slips using UV lithography. The typical diameter of the 
structures was 2 mm with 0.8 mm high walls. In closed 2D 
square shaped arenas we observed the roundabout motion of 
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FIG. 5: (a) Experimental snapshot of a circulating cell group with the velocities of cells in the corner of a 2 x 2 mm 2 square shaped 
microfabricated arena, (b) Simulations of model cells confined to a square arena show the emergence of circular motion over a wide range of 
model parameters. Supplemental video showing the circular motion of cells in the whole arena as well as in simulations are available online 



large cell groups. Fig. [5] shows an experimental snapshot of 
the circular motion with the instantaneous velocities of cells 
in the corner of a square shaped arena. A spectacular sustained 
whirling motion of the cells can be seen in the corresponding 
supplemental video B17I1 . 

We were also able to reproduce the above effects of bound- 
ary conditions on cell motion in simulations similar to those 
described in the previous section. Placing model cells in a 
square box with repulsive walls implemented through adding 
a repulsive force with an exponential falloff and finite a cutoff 
to the right side of equation (01 corresponding to each of the 
four walls 

Tp (j \_„ v / -Fwall cx P(~^rt) > d iw < Rq 

r w (ai W ) — n w x< "o ; 

L U , Kq < di w 

(6) 

where d- lvJ is the distance between cell i £ {1, N} and any 
of the w € {1,4} four walls with unit normal vector n w and 
-Pwaii = 50. In simulations cells circulated in an organized 
fashion under a wide range of noise and density values as 
well as for different system sizes. The online supplemental 
material II 1711 contains a simulation videos showing the emer- 
gence of organized circular motion of model cells confined to 
a square box. 

V. DISCUSSION 

In summary, we have presented evidence that purely short- 
range forces and simple experimentally motivated dynamics 



can be equivalent to an effective alignment term. Drawing an 
analogy between our experimental and model results, we im- 
ply that the emergence of collective motion among keratocytes 
is an example of a continuous kinetic phase transition. Our re- 
sults are also relevant in the broader context of recent work by 
Gregorie and Chate 01911 questioning the continuous nature of 
the transition in the angular noise model. Our present numer- 
ical results (exponents being similar to the ones determined 
for the original Vicsek et al. model H] ) support the view ll20ll 
that, in contrast with Ref. [19], angular noise models define a 
universality class with a corresponding continuous phase tran- 
sition in the ordering of the velocities. We expect that our and 
similar experiments as well as the quantitative model of the 
observations we have provided for the collective motion of 
tissue cells will lead to a better understanding of such vital 
phenomena as wound healing or embryogenesis. 
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